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j  Abstract 

\  This  paper  presents  approximations  for  the  rms  error  of  the  maximum 
likelihood  estimator  of  the  direction  of  a  plane  wave  incident  on  a 
random  array.  The  sensor  locations  are  assumed  to  be  realizations  of 
independent,  identically  distributed  random  vectors.  The  second  part  of 
the  paper  presents  an  asymptotically  unbiased  estimator  of  the  noise 


wavenumber  spectrum  from  random  array  data,  J, 
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Introduction 


Sonobuoy  fields  are  used  to  detect  submarines.  Thorn,  Booth,  and 
Lockwood*  have  proposed  that  the  signals  from  randomly  deployed  sonobuoys 
be  coherently  combined  to  make  acoustic  measurements.  They  present  the 
expected  value  and  variance  of  the  pattern  function,  and  the  distribution 
of  the  directivity  index  of  a  three-dimensional  random  array.  In  their 
model,  the  sensor  locations  are  observed  realizations  of  random  variables 
that  may  be  correlated  and  have  different  distributions.  They  define  an 
array  to  be  totally  random  if  the  sensor  locations  are  realizations  of 
independent,  identically  distributed  random  variables.  Several 
stochastic  properties  of  the  sidelobe  pattern  of  a  totally  random  array 
are  given  by  Steinberg. ^ 

The  ratio  of  the  peak  sidelobe  to  the  main  lobe  and  the  directivity 
index  of  an  array  system  are  measures  of  its  ability  to  perform  its 
tasks.  The  generic  signal  processing  tasks  of  an  array  system  are: 

1)  detecting  and  estimating  parameters  of  coherent  wave  signals  that 
impinge  on  the  array;  2)  resolving  multiple  wave  signals;  3)  estimating 
range,  bearing,  or  velocity  of  a  source  that  generates  the  detected 
signal;  and  4)  estimating  the  frequency-wavenumber  spectrum  of  the 
ambient  noise  field.  This  description  of  system  tasks  emphasizes  the 
statistical  nature  of  the  problem  of  measuring  performance,  especially 
for  random  arrays. 
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This  paper  presents  approximations  for  the  mean-square  error  of  the 
maximum  likelihood  estimator  of  the  bearing  of  a  plane  wave  impinging  on 
a  random  array  from  a  distant  source.  The  second  part  deals  with 
estimating  the  ambient  noise's  wavenumber  spectrum. 

1.  Random  Planar  Arrays 

Consider  a  planar  array  of  M  sensors  where  the  sensor  locations 
{(xk,yk)}  are  realizations  of  independent,  identically  distributed  random 
variables  {(Xk,Yk)}.  Assume  for  simplicity  that  the  signal  is  a  single 
frequency  plane  wave  plus  stationary,  zero  mean,  Gaussian  noise.  Let  0O 
denote  the  wave's  direction  of  arrival  with  respect  to  the  x  axis.  This 
angle  is  the  source  bearing  if  the  medium  is  horizontally  homogeneous. 

Let  u>0,  X0,  and  A  denote  the  wave's  frequency,  wavelength,  and  complex 
amplitude,  respectively.  The  signal  at  the  kth  sensor  is 

s(t,xk,yk)  -  AexpIKwQt-K^-Xyyjj)]  +e(t,xk,yk),  (1) 

where  kx  **  (2ir/X0)cos90  and  <y  m  (2iT/Xo)sin0o  are  the  x  and  y  components 
of  is  the  wavenumber,  and  e(t,xk,yk)  is  a  realization  of  the  noise 
field. 

The  correspondence  between  beamforming  and  frequency-wavenumber 

processing,  and  an  approximation  to  the  maximum  likelihood  (ML)  estimator 

n  M  _  9 

of  0O  have  been  presented  in  a  previous  paper. J  if  pEk«i(xk-x)z  and 
M  _2 

pZk-l'Yk-y)*  are  large,  where  p  is  the  power  signal-to-noise  ratio  in  a 

—  |  M  A 

narrow  band  about  uQ  and  x  ■  "  AEk-lxk»  LevinH  shows  that  the  root 


mean-square  errors  of  the  ML  estimators  of  kx  and  Ky  are  approximately 


-  M  _  "1/2 

rmse  kx  =  [2p  J  (x^-x)2] 
k-1 


M  _  -1/2 

rmse  <y  *  [2p  X  (yk~y)2] 
k-1 


(2) 


*  M  _  _ 

Moreover,  the  covariance  is  E(kx  -  xx)(iCy  -  <y)  “  (2pEjt»^(xjt-x) (y^-y)  )" 

These  expected  values  are  conditional  on  a  realized  array  geometry,  i.e 

they  are  ex-post  the  deployment  of  the  array. 

To  approximate  these  errors,  assume  that  M  Is  large.  Since  the 

sensors  must  lie  in  some  closed  and  bounded  set,  the  random  variables 

(Xk,Yk^  are  bounded.  Thus  the  central  limit  theorem  implies  that 

M  M 

M-^k-ltKk-)2  -  ox2  +  Op(M“1/2)  and  M-1!*.^-?)2  -  oy2  +  0p(M-l/2)t 

where  ox2  and  Oy2  are  the  variances  of  Xk  and  respectively,  and 

Op(M“1/2)  means  that  for  any  e>0,  there  is  a  Be>0  such  that  the  error  i 

bounded  by  B€M"l/2  with  probability  1-e .  Thus  the  rms  errors 

A  A 

of  kx  and  Ky  are  approximately 


rmse  kx  *  (2pM)"l/20“l 

x 

(3) 

A 

rmse  <y  *  (2pM)"l /2a-l 

y 


for  large  M.  The  estimators  are  approximately  uncorrelated  if  the 
coordinate  system  is  rotated  to  make  the  covariance  oTO  »  0  after 
rotation. 
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The  maximum  likelihood  estimator  of  the  bearing  Is 

AAA  A  A 

0O  ■  tan~^(<y/<x)  radians.  The  linear  approximation  of  tan~^(Ky/tcx)  - 
tan~^(Ky/fx)  Is 


-2  -1  * 
(l-hcv2icx  )~1[kv 


-2 


y  y' 


y*x 


(kx-*x)]. 


(4) 


Since  k  and  tc  are  approximately  uncorrelated  If  =  0,  It  follows 

a  y  xy 

from  (3)  and  (4)  that  when  pMox2  and  pM oy2  are  large, 

A 

E(0O-6O)2  =*  (X0/2it  )2(2pM)“*(a-2Bin20o  +  o~^cos^0o).  (5) 

x  y 

Thus  If  ox  ■  Oy  =  a,  then  from  (5) 

A  1/2 

rmse0o  “  X0(2pM)“  (2ira)_1rads.  (6) 

For  example,  let  a/XQ  «■  12,  M  -  90,  and  p  »  l/4(-6  dB) .  Then  from  (6), 
rmse0o  ■  0.11°  (1.98  x  10“-*  rads).  If  o/XQ  ”  100,  M  ■  40,  and 

A 

p  -  -10  dB,  then  rmse0o  *  0.03°. 

Now  suppose  that  and  are  independent  uniform  variates  whose 
range  Is  (0,L),  i.e.  the  sensors  are  uniformly  distributed  on  the  square 

A 

{0<x<L,  0<y<L}.  Then  a2  »  L2/12.  Let  us  compare  the  rmse0o  of 

O 

this  random  array  with  that  of  the  square  lattice  array  whose  M-N* 
sensors  are  at  the  points  { (jd,id) : j ,1-1 , . . . ,N) .  If  the  length  of  the 
square's  sides  Is  L,  then  the  sensor  spacing  is  d  ■  L/(N-1). 

From  (2),  (4),  and  (5),  we  only  have  to  compare  M-*£(xk_x)2  » 
M-^Cyk-y)2  with  o2. 


Since 
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M  N  _ 

M"1  I  Uk-*>2  -  M_1d2N  l  (j-j)2  (7) 

k-1  j-1 

,  L2  N+l 

-  d2(N-l)(N+l)/12  -  - 

12  N-l 

«  L2/ 12  -  a2, 

expression  (6)  holds  for  the  square  lattice  array.  The  approximate 
raise  of  the  maximum  likelihood  bearing  estimator  for  a  uniform 

A 

random  array  on  a  square  is  equal  to  the  approximate  rmse60  for  a 
uniformly  spaced  lattice  array  on  the  same  square. 

2.  Three-Dimensional  Random  Arrays 

For  a  given  coordinate  system,  let  Xk  =»  (xk»yk*zk)'  denote  the 
vector  location  of  the  kth  sensor  in  a  three-dimensional  array.  Let  0O 
denote  the  azimuth  angle  of  propagation  with  respect  to  the  x  axis,  and 
let  aQ  denote  the  elevation  angle  with  respect  to  the  z  axis.  Thus  the 
signal  at  the  kth  sensor  is 

s(t,Xk)  “  Aexp[i(u>0t  -  <.'xk))  +  e(t,Xk>» 

where  jc '  -  (xx,<y,<z)  is  the  vector  of  wavenumber  components 
kx  *  (2x/Xo)cos0o,  <y  -  (2ir/Xo)8in0o,  and  <z  -  (2n/X0)coso0. 

The  correspondence  between  beamforming  and  frequency-wavenumber 
processing  holds  in  three  dimensions.  The  ML  estimators  of  the 
wavenumber  components  are  the  kx,  ,  and  <z  that  maximize 


T 
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N  M 

I  l  l 


j-1  k-1 


s(bj  »xk*yk»zk^exP^^ii,2Ek  "  uot j  ) J 


2 


(8) 


where  N  is  the  number  of  simultaneous  discrete-time  observations  of  the  M 

A  A 

channels.-*  The  rms  errors  of  <%  and  <y  are  approximated  by  (2),  and 
*  M  _  y 

rmseKz  a  l2pSjCJBj(zjt-z)^]-1/2  .  once  again,  the  ML  estimator  of  the  source 

A  A  A 

bearing  is  0Q  -  tan~^(<y/Kx) ,  and  thus  (5)  holds  for  a  totally  random 
three-dimensional  array  of  M  sensors. 

3.  Estimating  the  Wavenumber  Spectrum 

Consider  the  problem  of  estimating  the  frequency-wavenumber  spectrum 
of  the  ambient,  zero  mean,  Gaussian  noise  field  around  a  random  array. 
Since  an  n-dimensional  array  is  not  much  harder  to  analyze  than  a  linear 
array,  let  x^  ■  , . . . *  denote  the  vector  position  of  the  kth 

sensor  with  respect  to  a  fixed  coordinate  system.  Assume  that  the  x^  are 
realizations  of  independent  random  vectors  {Xj^  *  (Xk^ , . . . ,X^n) '}  that 
have  a  common  continuous  multivariate  density  f (x) .  Rotate  the 
coordinate  system  so  that  the  covariance  matrix  of  X^  is  diagonal,  and 
for  simplicity  let  Oj2  «  ...  «  an2  =.  a2^  a2  tjje  variance  Qf  each 

X^£  after  rotation. 

Let  e(t,x)  be  the  noise  at  point  x  at  time  t.  If  the  noise  field  is 
stationary  in  t  and  jc,  the  covariance  function  ce(t,£)  =  Ee (t+T »x+jOe (t ,x) 
is  independent  of  t  and  x.  The  frequency-wavenumber  spectrum  is  defined 
as 


S€(w,jO  -  /c€(T,2)exp[i(£'£  -  u)T)ld^_, 


(9) 


1 
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assuming  that  ce  Is  absolutely  integrable.  The  power  spectrum  of  the 
noise  is  Se(u>,0). 

Assuming  that  the  channels  are  sampled  at  times  tj*jA  for 

j«*0,...,N-l,  define  the  discrete  Fourier  transform 
N-l 

**  e(jA)exp(-iwJA) :  k  -  1,...,M).  If  Sc(u>,0)  is  bandllmlted 

at  u/A,  then  N“^E|c(xjc)  | ^  «  A“*Se(<»>*0)  for  large  N.^  Let  us  work  with  the 
e(xfc)  to  obtain  an  estimator  of  Se(jc,o>)  for  a  given  10,  which  will  be 
denoted  SeOO  to  simplify  notation.  The  properties  of  the  estimator 
depend  on  the  following  theorem. 


Theorem.  Define  the  n-dimensional  Fourier  transform,2 
UOO  ■  1^*1  eCx^Jexpdjc'x^) .  Assume  that  D(o)  -  ff2(x)dx  -  0(o“n)  and 
whence  *  0,  |$0e)|  <  co-11  for  some  constant  c,  where  $(<)  -  Eexp(ijc ’X^)  is 
the  characteristic  function  of  These  assumptions  hold  for  the 

multivariate  normal  and  uniform  densities.  Then 

lim  (DM2)_1e|u(k)|2  -  Se(*0, 

M,o-m»  “ 

and  U(kj)  and  U(xy)  are  asymptotically  uncorrelated  for  jtj  *  <9. 

M 

Proof:  The  array  transfer  function  is  R(jc)  «  E^.^expCi^'x^) .  For  large 
H,  M~*R(k)  ■  $00  +  0p(M“1/2)  by  the  central  limit  theorem.  Thus 

(DM2)-1R(<1)R*(k2)  -  D_^(kj)^*(£2>  +  0p(M-l/2)  (10) 

(star  denotes  complex  conjugate)  since  D“ 1 1 ^ ( k_)  |  -  0(1)  in  the  cross 
product  by  the  above  assumptions.  Thus 


'& 

* 
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lim  (DM2)-1(2it  )~nf  |R(<)  |  2  dK  -  D— ^ ( 2n  )“n/  j <J> ( *c )  |  2d< 

M+oo  ~  ~ 

=  D*1Jf2(x)dx  -  1.  (11) 

From  (10),  lim^*,  (DM2)-1  |R(0)  j2  =  D-1 14>(0)  | 2  -  D-1  =  0(an).  Thus  (11) 
implies  that  as  M  and  a-*-*,  (DM2)-^ | R(£) | 2+6(£) ,  a  Dirac  delta  function.  If 

£l*£2» 

(DM2)_^R(k^1)R*(£2)  -  o (o“n)  +  0p(M-l/2).  (!2) 

These  limit  results  are  used  as  follows: 

M  M 

E[U(kj)U*(£2)]  “I  l  ce(*l  “  xk)exp[i(£i,Xj  “  £2'Sc^  (13) 

3=1  k=l 

M  M 

=  (2tt  )“n  l  l  /sc(v)exp[-iv(x.j  -  x^)  JexpUOci'Xj  -  <2'^) 
j=l  k=l 

from  the  inverse  of  (9).  Gathering  terms, 

E[U(kj)U*(<2)]  =  ( 2tt ~  v)R*(<2  “  ^se^)L^d3L*  04) 

Thus  from  the  above  limits  and  (14),  limM>a+oo(DM2)''1E|uO<)  |2  = 

(  2ir  )~nJ 6 (x-v)  Se  (v)d_v  =  Se  (<_) . 

If  *  jc 2 ,  then  lim^  rlA.CB(DM2)~^E[U(ic ^ )U*(<?)  1  =  0  from  (12).  Thus 
U(£i)  and  U(<?)  are  asymptotically  uncorrelated.  For  finite  M<<o2n,  the 
correlation  is  Op(M”l/2). 

This  theorem  provides  a  basis  for  estimating  Se(j<).  One  method  is  to 
divide  the  (time)  sample  into  J  segments  of  successive  observations, 

Nj  ■  N/J,  and  compute  U(£)  for  each  segment.  These  UjOO's  will  be 
approximately  uncorrelated  if  Nj  is  large.  Thus  from  the  theorem, 


9 

Uj (O  have  a  complex  Gaussian  distribution  for  each  j  (the  noise  Is 

Gaussian),  2(DM2)“1  |  Uj  (<_)  |  ^/Se(0  is  approximately  chl-squared  with  two 

degrees  of  freedom  and  thus  the  variance  of  SeOO  is  approximately 

J_1S2(k). 
e  — 

4.  A  Planar  Array  Example 

Continuing  with  the  vector  notation,  suppose  that  the  sensors  are 
uniformly  distributed  on  the  square  {-L/2  <  x^  <  L/2,  -  L/2  <  X2  <  L/2} . 
Thus  f(x)  ■  1/L2  for  x  In  the  square,  *  0£2  =*  a2  *  L2/12,  and 

D  *  /f2(x)dx  *»  L“2.  The  assumptions  for  the  theorem  hold  since  D  =  0(o“2) 
and  $(jc)  *  4(x1X2L2)“lsin(KiL/2)8in(K2h/2)  =  0 (a-2).  Thus 
(L/M)2E|U(0 |2  -  SeCO  for  large  M  and  L  in  this  example.  The  estimator 
of  Se(0  is  then  (L/M)2J“1Ej1.1|Uj(0  [2  using  the  time  segmentation  method. 


*Thi8  work  was  supported  by  the  Office  of  Naval  Research  (Statistics  and 
Probability  Program)  under  contract. 
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on  the  grid,  then  the  FFT  algorithm  can  be  used  to  compute  U(k). 


* 

J 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  gnltrxQ 

I  REPORT  DOCUMENTATION  PAGE 


/\/  j  '  /  t/^ 


REPORT  NUMBER 


Technical  Report  22 


a 

dp 


2.  GOVT  ACCESSION  NO. 


L«.  TITLE  r«nd  Submit) 


READ  INSTRUCTIONS 

_ BEFORE  COMPLETING  FORM 

"1  recib*e'rt,s  catalog  NUMBER 


ORT  ft  PERIOD  COVERED 


Estimating  Signal  and  Noise  using  a  Random  Array 


autmo«c«; 


[Melvin  J J  Hinich 


: Tecnnica  1  /epflrt 


6.  PERFORMING  ORG.  REPORT  NUMBER 


E  jCORT^ACT-.OR  GAANT  TTUMBERr*^ 


j  Hini< 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  7  .  " 7 

Virginia  Ta«h  I  ;  l  ,  f 

Department  of  Economics  v/ 

Blacksburg,  VA  24061 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research  Code  426  f  t 

Statistics  and  Probability  Program  / 

Arlington,  VA  2221  7  _ _ _^2 

TT  MONITORING  AGENCY  NAME  a  ADDRESS#!/  dllltfnt  from  Controlling  Olllct) 


:  i  j 


MW14-75-C-9494  v  J 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  «  WORK  UNIT  NUMBERS 


NR  042-315 

r 2 ■  REPORT  DATE  - - 1 

I  7jul**#81  / 

•13  ?  MUM0ER  OF  — — — 

- »  11 _ 

15.  SECURITY  CLASS,  (of  thia  rmport) 

Unclassified 


15a.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thi*  R#pofT) 

Distribution  of  this  document  is  unlimited 


[  17.  DISTRIBUTION  STATEMENT  (oi  tha  abstract  antarad  In  Block  20.  II  dtlferant  from  Rmport) 


18.  SUPPLEMENTARY  NOTES 


19-  KEY  WORDS  (Contlnua  on  ravaraa  a  I  da  It  nacmaamry  and  Idantlty  by  block  numbar) 

Random  Array,  Array  Processing,  Wavenumber  Spectrum,  Array  Response, 
Bearing  Estimation 


1 20.  ABSTRACT(Contlnuaonravaraaaldaltn0caaaary  and  Idantlty  by  block  numbar) 


DO  ,  :°r73  1473  EDITION  OF  *  NOV  69  IS  OBSOLETE 

S/N  0102-  LF-014-  6601 


HAQJilO 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  flFhftn  Dftf*  BntftrftO 


